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In this Letter we describe a novel class of dynamical excitations - accelerating oscillatory fronts 
in a new genre of nonlinear sonic vacua with strongly non-local effects. Indeed, it is surprising that 
such models naturally arise in dynamics of common and popular lattices. In this study, we address 
a chain of particles oscillating in the plane and coupled by linear springs, with fixed ends. When 
one end of this system is harmonically excited in the transverse direction, one observes accelerated 
propagation of the excitation front, accompanied by an almost monochromatic oscillatory tail. The 
front propagation obeys the scaling law I ~ The frequency of the oscillatory tail remains 

constant, and the wavelength scales as A ~ t^^^. These scaling laws result from the nonlocal effects; 
we derive them analytically (including the scaling coefficients) from a continuum approximation. 
Moreover, a certain threshold excitation amplitude is required in order to initiate the front propaga¬ 
tion. The initiation threshold is rationalized on the basis of a simplified discrete model. This model 
is further reduced to a new completely integrable nonlinear system. The Letter introduces a new 
and yet unexplored class of nonlinear sonic vacua and explores the effects of strong non-locality on 
the initiation and propagation of oscillating fronts in these media. Given their simplicity, nonlinear 
sonic vacua of the type considered herein should be common in periodic lattices. 


Nonlocal nonlinearities naturally appear in classical 
models describing the nonlinear motion of beams with 
restrained ends [1, 2]. For example, the transverse dy¬ 
namics of a beam with fixed ends, bending in a direction 
of one of the main axes of its cross-section, is described 
by the following partial integro-differential equation: 
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Here y{x,t) is the transverse displacement of the beam 
point; a; is a coordinate along its axis; p and E are the 
mass density and Young’s modulus, respectively, of the 
beam; L is the length of the undeformed beam; S and 
I denote area and moment of inertia of the beam cross- 
section and T is an applied axial tensile force in x direc¬ 
tion. The nonlocal term in the right-hand side of Eq.(l) 
appears due to a stretching of the midplane caused by 
transverse displacement of the beam with immobile ends. 
If the bending term is negligible, the axial tensile force is 
absent, and both ends are clamped at distance L; then 
Eq.(l) reduces to the following strongly nonlinear wave 
equation: 
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This equation describes the transverse oscillations of an 
elastic string with fixed ends and without pre-tension. 
Only recently, a similar approximation has been devel¬ 
oped for a discrete counterpart of such a string in the 
form of a linear chain of particles moving in the plane, 
with next-neighbor interactions through linear springs, 
and with fixed boundaries [3]. As shown in [3] and Sup¬ 
plemental Material [4] , if the transverse displacements of 
the chain of particles are not too large compared to the 
length of the free springs, the transverse motion of the 
chain of particles decouples from the longitudinal mo¬ 
tion, and is approximately described by the following set 
of ordinary differential equations: 
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In this system yn is the transverse displacement of the 
n-th particle, N - the overall number of the particles, 
yi = yN = 0. The mass of each particle and the spring 
constant are set to unity without loss of generality. It 
is easy to recognize that Eq. (2) approximates Eq.(3) in 
the long-wave limit. 

Systems (2) and (3) exemplify the important concept 
of sonic vacuum - linearization of both these systems 
yields zero sound velocity. Thus, these systems can be 
classified as nonlinearizable - one never can neglect the 
nonlinear terms. Similar properties of the sonic vacuum 
are well-known and widely studied in systems of granular 
particles without external pre-compression; the nonlin¬ 
earity appears there due to Hertzian contact [5-8]. Also, 
nonlinearizable systems of different structure were widely 
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explored as possible nonlinear energy sinks [9-12]. The 
essential nonlinearity of these sinks allows them to res¬ 
onate with primary oscillatory systems over a broad fre¬ 
quency range [10] and efficiently absorb energy. Such 
systems are investigated as possible engineering solutions 
for vibration mitigation and energy harvesting in a wide 
range of possible applications [10, 13, 14]. 

However, systems (2), (3), as compared to the granu¬ 
lar media without pre-compression or the nonlinear en¬ 
ergy sinks, possess an additional important property - 
the nonlocality. If one is interested only in modal oscil¬ 
lations of these systems, the problem becomes relatively 
easy if one notes that the integral term in Eq. (2) and 
the corresponding sum term in Eq. (3) depend only on 
time. Therefore, the spatial modal shapes of the strongly 
nonlinear sonic vacua (2) and (3) will be the same as for 
simple linear string or chain with fixed ends, and, respec¬ 
tively, will correspond to sinusoidal standing waves. The 
nonlinearity will reveal itself only in the time domain - 
the oscillations will be anharmonic, and their frequency 
will be proportional to the amplitude [3]. The fact that 
the integral term in Eqs. (1), (2) does not affect the 
spatial modal shapes, is well-known and widely used in 
approximate modal analysis of oscillating continuous sys¬ 
tems [15]. 

Beyond modal oscillations, the dynamics of sonic vacua 
similar to (2), (3) is almost unexplored due to the strong 
nonlinearity. Below we demonstrate that the nonlocal 
nonlinearity brings about some unexpected and unusual 
dynamic phenomena even in very simple settings. To il¬ 
lustrate that, the dynamics of a spring-and-mass chain 
with N particles in a plane is simulated. This system 
is detailed in the Supplemental Material (system (1), 
[4]), and in the limit of low energy its transverse oscil¬ 
lations are approximately described by Eq. (3). The 
non-dimensional Hamiltonian of the lattice is expressed 
as 

Af-l N-l 

n—2 n—1 

where Xn and denote the axial and transverse coordi¬ 
nates of the of the n-th particle, r„ = [(x„+i — -I- 

{Un+i — 2/ti)^]^^^ is the distance between particles n and 
n -I- 1. The right end of the chain is fixed, whereas the 
left end of the chain does not move in the x direction, 
and harmonically oscillates in y direction. The equations 
of motion are expressed as 

xi(t) = 0, yi{t) = Hsin(wt), 

Xn = -dH/dxn, yn = -dH/dyn, n = 2,...,N -1,{5) 

XNit)=N-l, 2/Ar(t)=0. 

and are numerically integrated with zero initial condi¬ 
tions by the velocity Verlet method [16] with the follow¬ 
ing parameters: N G [125,10000], excitation amplitudes 
A G [0.1,7] and frequencies uj G [0.01,0.5]. Small vis¬ 
cous friction is imposed on a few rightmost particles to 




FIG. 1: Evolution of local kinetic energy E = y^/2AE at 
different segments of time t — t/At in the chain with N = 625 
and (a) A = 0.5, uj = 0.05 (AE = 10"®, At = 492); (b) 
A = 1, UJ = 0.3 (AE =l,At= 164); (c) A = 0.15, uj = 0.05 
(AE = 10“^, At = 2048). 


improve the simulation accuracy and avoid numerical dis¬ 
crepancies. 

Typical results of the simulation are presented in Fig¬ 
ure 1, where we depict the ’’transverse component” of 
local kinetic energy Bk = y^l2 versus n at different time 
instances. One can observe the propagation of the exci¬ 
tation front, accompanied by an apparently monochro¬ 
matic oscillatory tail. Moreover, it is clear from Figure 1 
(a), that this front accelerates in the course of propaga¬ 
tion. It is interesting to note that the stationary fronts 
with oscillatory tails, but with constant velocity, are well- 
known in models of phase transitions in solid state and 
similar problems [17-20]. Moreover, it is obvious that in 
a linear wave equation with similar boundary excitation 
one would observe the monochromatic oscillatory front 
propagating with the sound velocity. It follows that the 
observed acceleration of the front should be attributed to 
the nonlocality of the nonlinear term in (3). In order to 
explain this finding analytically, we consider a simplified 
model of the oscillatory region in the chain and suppose 
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a monochromatic wave in the oscillatory tail after the 
front in continuum approximation. Only transverse os¬ 
cillations are taken into account (numerical justihcation 
of this assumption is presented in the Supplemental Ma¬ 
terial [4]). Then, the field of displacements in the oscil¬ 
latory zone is described as follows: 


y{x,t) 


As\n{ijjt — kx), 0<x<l{t), 
0, l{t) < X < L. 


Here l(t) is the instantaneous coordinate of the front, and 
k is the wavenumber. It is also assumed that the front 
propagation is slow enough compared to the frequency 
of transverse oscillations of the particles; i.e., dl{t)/dt 
wl{t). To establish complete correspondence between the 
continuum approximation (2) and the discrete model (3)- 
(5), one should set p = E = S = 1, L = N — 1. Then, 
by substituting (6) to (2) and balancing principal terms, 
the following equation is obtained: 

. 9 . / . N k‘^A^l{t)s\niujt — kx) 

sin(wt - kx) ^^ 

An additional condition can be obtained from the as¬ 
sumed stationary character of the front propagation. To 
this end, the phase velocity of the oscillatory tail should 
be equal to the front velocity [19] so that Vph = uj/k = 
dl{t)/dt. Combining this expression with Eq.(7), one ob¬ 
tains an explicit expression for the position of the accel¬ 
erating front: 


^ A2 ^ l^dt’ ~ 4{N-1) 
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l{t) = K = 
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Thus, the front indeed accelerates with velocity V ^ 
Prediction of Eq. (8) is completely supported by the nu¬ 
merical simulations, as is demonstrated in Figure 2 - for 
three different sets of parameters, with the curves de¬ 
picting the front position versus time ( shifted by In K) 
collapsed into a straight line with slope 4/3. So, the con¬ 
siderations presented above predict not only the correct 
scaling law for the front position, but also the scaling 
coefficient, which depends on specific set of parameters. 

The simulated system is discrete rather than continu¬ 
ous, and that is why not every set of parameters leads to 
formation of the accelerating front. This point is illus¬ 
trated in Fig. 1 (b): If the excitation amplitude is too 
small, the oscillations remain localized at the left end 
of the chain. This phenomenon cannot be explained in 
terms of the continuum model (2). To describe the front 
initiation above a certain excitation threshold, we resort 
to the discrete model of the chain. We will adopt a sim¬ 
plified approach and establish the minimal amplitude of 
oscillations of particle n = 1 that allows efficient excita¬ 
tion of particle n = 2 and, thus, substantial excitation of 
the chain and initiation of the wave front. Accordingly, 



FIG. 2: (Color online) Dependence of the position of the lead¬ 
ing edge of the front l{t) on time t for three different sets of 
parameters: N = 2500, A — 0.5, ta = 0.03; N = 5000, A — 
1.5, ta = 0.08 and N — 1000, A = 0.6, ca = 0.03 (curves 1, 
2 and 3 respectively). Dashed line corresponds to the slope 



FIG. 3: (Color online) Phase trajectories correspond¬ 
ing to Equation (12)for zero initial conditions and a = 
0.3,0.34,0.3575,0.37,0.42 (curves 1, 2, 3, 4 and 5 respec¬ 
tively) . 


we analyze (3) with yi = —Asina;t,?/„ = 0,n > 3 and 
zero initial conditions for y 2 {t)- This system is rescaled 
with ojt = T, ?/2 = aulr), A = aa, a = ujy/2(N — 1). 
Then one arrives at the following equation for variable 
m(t): 

Utt + (2u -I- asinr)[u^ -I- (it -I- asinr)^] = 0. (9) 

The primary frequency of the oscillatory front is expected 
to be close to the normalized value of unity. Therefore, 
the complex variable (^(r) exp(iT) = Urir) + iu{T) is in¬ 
troduced [21]. Supposing that variable if{T) varies slowly; 
balancing principal terms in Eq. (9), we arrive at the fol¬ 
lowing slow-flow equation: 

V>r + - ^*[2(2(p -b a)[\ip\^ + \^ + a\^) + (10) 

+ a)(i^^ + ((/? + q)^)] = 0. 
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the chain, similar to the regime demonstrated in Fig. 1 
(b). There exists an excitation threshold above which 
the phase trajectory is attracted to a region of relatively 
large R. So, the energy of oscillations is intensively irra¬ 
diated into the chain, and it seems natural to associate 
this regime with the formation of the oscillatory front. 
The threshold excitation corresponds to the phase tra¬ 
jectory, which passes through the saddle point in Fig. 3 
(curve 3). This yields the following evaluation for this 
threshold: 


a^r ~ 0.3574 


(13) 


FIG. 4: (Color online) Correlation between the frequency uj 
and the amplitude A for the wave front for chain lengths 
N = 125 and N = 250 (lines 1 and 2). The black dashed 
lines represent linear fitting of the markers 1 and 2. The 
region above a corresponding dashed line corresponds to a 
propagating front regime, and no front is initiated for pa¬ 
rameters below the dashed line. The linear fitting rela¬ 
tions are given hy A = 0.274a; y^2(A'" — 1) for N = 125 and 
A = 0.287a;Y^2(iV - 1) for N = 250. 

Though far from obvious, this slow-flow equation is com¬ 
pletely integrable. The integral of motion is expressed 
as: 

C = + \<P + + \<p'^ + {^ + a)'^\ ]• ( 11 ) 

To see this, it is sufficient to note that Eq. (11) is 
equivalent to Lp-r = — Since C is real, one im¬ 
mediately obtains consequently, ^ = 

= 0- Then we split the slow variable 
into polar components ip = Rexp(iS). Initial condition 
:p(0) = 0 corresponds to C = —3a^/8. Therefore, the 
averaged phase trajectories of the particle n = 2 for dif¬ 
ferent values of the external excitation are expressed by 
the following family of implicit equations: 

R^{l — ^R'^ — a'^) — 3Ra{R^+^a^)cosS—2{Racos6)'^ = 0 . 

( 12 ) 

The family of solutions of Equation (12) for various 
values of a is presented in Fig. 3. One can see that for 
small values of a the phase trajectory stays in the re¬ 
gion of small values of R. This regime corresponds to 
localization of forced oscillations near the excited end of 


So, the boundary for formation o f the oscil latory front is 
described by the line A = acr<^^j2{N — 1). This predic¬ 
tion is verified in Figure 4. Approximate linear depen¬ 
dence of A on w is observed, but the coefficient of these 
lines turns out to be somewhat overestimated. 

This discrepancy is obviously related to the number 
of simplifying assumptions adopted in our analysis. Be¬ 
sides, it is apparent from Fig. 1 (b), that localized oscilla¬ 
tions at the left end of the chain exhibit chaotic dynam¬ 
ics. Then, it is possible to expect that the appropriate 
conditions for the front initiation may be formed also if 
the excitation is far below the threshold - just due to 
fluctuations. An example of such behavior -- the front 
formation is observed after a certain time delay - is pre¬ 
sented in Fig. 1 (c). 

To conclude, we revealed a new type of excitations 
in a lattice representing a nonlinear sonic vacuum with 
strong nonlocal dynamical interactions (despite only 
next-neighbor physical coupling). These excitations are 
accelerating fronts with oscillatory tails. The fronts ac¬ 
celerate according to the scaling law I ~ due to non¬ 
local dynamical interactions. The tails have constant fre¬ 
quency, but their wavelength is not constant - it scales 
with time as A = 27r//c = 2'Klt/w ^ Such fronts re¬ 
veal themselves in most well-known and popular models, 
such as the suspended string without pre-tension and the 
chain of linear springs and masses with fixed ends. Simple 
analytic considerations allow derivation of all main pa¬ 
rameters of the front, including the scaling characteristics 
and the excitation threshold. Due to the fixed bound¬ 
ary conditions, such fronts can exist only as transient 
regimes. At the same time, due to extreme simplicity 
and popularity of the involved models, one can expect to 
see such accelerating fronts in many physical settings. 
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